Evolution of symmetry and structure of the gap in Fe-based superconductors with 

doping and interactions 
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We present a detailed study of the symmetry and structure of the pairing gap in Fe-based super- 
conductors (FeSC). We treat FeSC as quasi-2D, decompose the pairing interaction in the XY plane 
in s-wave and d-wave channels into contributions from scattering between different Fermi surfaces 
and analyze how each scattering evolves with doping and input parameters. We verify that each 
interaction is well approximated by the lowest angular harmonics. We use this simplification to 
analyze the interplay between the interaction with and without spin-fluctuation components, the 
origin of the attraction in the s and d x 2_ y z channels, the competition between them, the angular 
dependence of the s gaps along the electron Fermi surface, the conditions under which s gap de- 
velops nodes, and the origin of superconductivity in heavily electron- or hole-doped systems, when 
only Fermi surfaces of one type are present. We also discuss the relation between RPA and RG 
approaches for FeSC. 
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I. INTRODUCTION 

The symmetry and the structure of the superconduct- 
ing gap in Fe-based superconductors (FeSC), and their 
evolution and possible change with doping are currently 
subjects of intensive debates in the condensed matter 
community. 

The vast majority of researchers believe that supercon- 
ductivity in FeSC is of electronic origin and results from 
the screened Coulomb interaction, enhanced at partic- 
ular momenta due to strong magnetic fluctuations 1-13 
or orbital fluctuations 14-17 . For systems with a single 
Fermi surface sheet, such interaction cannot lead to a 
simple s-wave superconductivity, but it can give rise to a 
superconductivity with non-s-wave symmetry - p-wave 
for strong ferromagnetic spin fluctuations and d-wave 
for strong antifcrromagnctic spin fluctuations. In FeSC, 
however, the electronic structure is more complex. The 
low-energy states are formed by the hybridization of all 
five Fe-d orbitals what in the band basis not only gives 
rise to multiple sheets of the Fermi surface (FS) [weakly 
doped FeSC contain two electron FSs and either two or 
three hole FSs], but also leads to a complex mixing of 
contributions from intra- and inter-orbital terms in the 
interactions between low-energy fermions. In this situ- 
ation, in addition to the potentially non-s-wave super- 
conductivity, FeSC may also develop superconductivity 
with an s-wave symmetry of the gap even for the repul- 
sive electron-electron interaction. For parameters used 
for orbital interactions in most studies of FeSC, the s- 
wave gap, averaged over the FSs, changes sign between 
different FS sheets (it is commonly called the s gap), 
and the superconducting state competes with the spin- 
density-wave (SDW) state. But for other parameters, 



a conventional superconducting state, with a sign pre- 
serving s ++ gap, becomes possible, and such a state will 
compete with the charge-density-wave (CDW) state 11 . 

The existing theoretical approaches to pairing in FeSC 
can be broadly divided into two categories. One assumes 
that fermions at energies smaller than a fraction of the 
bandwidth can be treated as itinerant, with a moder- 
ate self-energy 18 , although strong coupling effects, such 
as interaction-driven renormalization of the whole band- 
width, have to be incorporated. 19-21 In the itinerant ap- 
proach, the pairing is often treated in a BCS/Eliashberg 
formalism, with the interaction taken as a combination of 
direct electron-electron interaction and effective interac- 
tion mediated by collective bosonic excitations. Another 
approach assumes that the system is not far from the 
Mott regime, and that the pairing should be affected by 
the tendency towards a Mott insulator 22-24 . 

This work falls into the first category. Already within 
this category, there are several, seemingly different ap- 
proaches to the pairing: the RPA-typc spin-fluctuation 
approach (RPA/SF) 1-3 ' 5 ' 6 ' 25-28 , the functional RG ap- 
proach (fRG) 7-9, 29,30 , and the analytic (logarithmic) 
RG approach based on a "minimal" model for the 
pnictides 10-12 , in which the interaction in each pairing 
channel is restricted to the leading angular harmonics in 
each pairing channel (the leading angular harmonics ap- 
proximation, or LAHA). The very positive fact for the 
itinerant approach as a whole is that, so far, the results 
of all these different approaches agree on the pairing sym- 
metry and the gap structure in hole-doped and electron- 
doped FeSC. Namely, all three approaches predict that 
the leading pairing instability at small/moderate dop- 
ings is in the s-wave channel, and the gap, averaged over 
the individual FS sheets, changes sign between hole and 
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electron sheets (an s ± gap). The s ± gap generally varies 
along each of the FSs. If the FeSC are treated as 2D 
systems (i.e., if the variation of the interaction along k z 
is neglected) , the variation of the gap is stronger on the 
electron FSs, and can be large enough to create nodes. 
All three approaches also predict that the gap with the 
nodes is more likely in either undoped or electron-doped 
FeSC, while in hole-doped FeSC a no- nodal state is more 
likely. 

The RPA and fRG formalisms have been also 
applied 26,30 ' 31 to study superconductivity in recently dis- 
covered heavily electron-doped KFe2Se 2 , where only elec- 
tron FSs remain, according to recent ARPES studies 32,33 . 
The results of RPA/SF and fRG approaches for this lim- 
iting case are again in agreement - both predict that 
the gap should now have d x 2_ y 2 symmetry, i.e., it should 
change sign between the two electron FSs. Other ap- 
proaches, however, found a conventional s-wave super- 
conductivity in this limit, by one reason 34,35 or the 
other 15 . A more complex order parameter has also been 
suggested. 36 

The d— wave gap symmetry was also predicted by 
fRG 29 for heavily hole doped KFe2As2, where only hole 
FSs are present. The RPA/SF analysis for this mate- 
rial 27 found attraction of comparable strength in both 
s-wave and d-wave channels. 

The goal of this paper is to understand in more detail 
the evolution of the s ± gap with hole and electron dop- 
ing and the interplay between the s^ and d x 2_ y 2 pairing 
in FeSC. The idea is to "decompose" the full pairing in- 
teraction into contributions from scattering processes be- 
tween different FSs and check how each process evolves 
with doping and input parameters. Such an analysis has 
been performed within the RPA/SF formalism in Rcf. 1 
but only for an s-wave interaction and only for a fixed 
doping and a particular set of input parameters. 

In this work, we combine the 5-band RPA/SF and 
LAHA approaches. Specifically, we take the full set of 
interactions TijQc, k') from the 5-band RPA/SF calcu- 
lation as input (i, j label different FSs, k and k' are 
momenta along these FSs) and show that, for all cases 
that we studied, different interactions r^^k, k') are well 
approximated by only two angular harmonics in s-wave 
and d x 2_ y 2 channels (we choose one of the momenta such 
that only these two pairing channels contribute). This 
leaves us with a finite number of interaction parameters 
in each of the two pairing channels (the number of inde- 
pendent parameters is 4 or 5, equal to the total number 
of FSs). We then solve the pairing problem, obtain the 
cigcnfunctions (which determine the gap structure) and 
eigenvalues (which are the dimensionless couplings) and 
analyze how both eigcnfunctions and eigenvalues evolve 
with the parameters. In this work, we neglect potential 
new physics associated with 3D effects and treat FeSC as 
quasi-2D systems, i.e., neglect the k z dependence of the 
quasiparticle dispersion and of the interactions. 

The key goal of our work is to understand whether su- 
perconductivity in FeSCs is governed by a single underly- 



ing pairing mechanism for all hole and electron dopings, 
despite that the pairing symmetry and the gap structure 
may change, and whether the entire variety of pairing 
states can be adequately described within the effective 
low-energy model with small numbers of input parame- 
ters. 

The specific set of issues that we address are: 

• What is the origin of the strong angular dependence 
of the s^ gap along the electron FSs? The angular 
dependence of the effective interaction Tij(k, k') 
due to the change in orbital character of the 
states on the FS, the competition of the scattering 
between hole-electron and electron-electron-sheets 
and the local Coulomb repulsion are candidates 
that can give rise to a strong anisotropy of the gap. 

• Are the angular dependencies of all interactions rel- 
evant for the gap structure, or can some interac- 
tions be safely approximated as angle-independent? 
In principle, the angular dependencies of both the 
electron-hole and electron-electron interactions can 
affect the structure of s gap. In the LAHA ap- 
proach, we can vary the angular dependence of each 
interaction by hand and explore how it affects the 
gap structure. 

• Why do the s solutions obtained within the 
RPA/SF and fRG approaches have nodes for sys- 
tems with two hole and two electron FSs and no 
nodes for systems with three hole and two electron 
FSs? Is this behavior generic, or just a trend mean- 
ing that in both cases the s ± gap is either nodal 
or non-nodal, depending on the input parameters? 
This issue is difficult to address in both RPA/SF 
and fRG approaches as these are numerical meth- 
ods which require certain input parameters and ex- 
ploring parameter space is computationally expen- 
sive. But it can be addressed within LAHA as one 
can continuously change any of the input parame- 
ters. 

• What causes the pairing when only electron FSs are 
present? Possibilities include pairing driven by the 
angular dependences of the interactions; s— wave 
pairing caused by virtual scattering to gapped hole 
states -"incipient" s-wave (s ± gap but without hole 
FSs), a d— wave pairing (a plus- minus gap on elec- 
tron FSs) due to magnetically enhanced repulsive 
interaction at momentum transfer (it, tt) between 
the electron FSs, or an s— wave pairing if the inter- 
action at (tt, tt) is strong and attractive. 

• What causes the pairing at large hole doping, when 
only hole pockets are present? It is possible that 
the scattering between hole FSs with wave vector 
(tt, tt) favors s pairing while the alternative is that 
Coulomb avoidance within each FS pocket makes 
nodal d-wave the preferable symmetry. 
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• How is the structure of the pairing interaction af- 
fected when the spin-fluctuation component is added 
to the direct fermion-fermion interaction'? Inclu- 
sion of the SF contribution can affect the relative 
magnitudes of the interactions r^k, k') between 
different FS sheets, and in principle can change its 
angular dependence. 

This last issue is relevant for understanding the com- 
parison between RPA/SF and RG approaches, which we 
pause to discuss in some detail. These two approaches 
differ in the assumption of what are the relevant energy 
scales for magnetism and superconductivity. In the RG 
approach it is assumed that magnetism and supercon- 
ductivity are produced by the same low-energy fermions 
and have to be treated on equal footing, starting from 
a model with only density-density and exchange inter- 
action between fermions. Then the SF component of 
the interaction develops together with the pairing ver- 
tex (one builds up stcp-by-stcp "parquet" rcnormaliza- 
tions simultaneously in the particle-hole and particle- 
particle channels) 7,11,12 . In the RPA/SF approach 2,5,6,25 
the assumption is that magnetism is superior to super- 
conductivity and comes from fermions at energies com- 
parable to the bandwidth, which are above the upper 
edge for the RG treatment. If so, the SF component 
should be included into the bare interaction, which does 
not need to be further renormalizcd in the particle-hole 
channel. In this situation, only renormalizations in the 
particle-particle channel (which always come from low- 
energy fermions) remain relevant, and the RG treatment 
becomes equivalent to the BCS theory. There is no good 
justification to select a particular set of diagrams for the 
renormalization of the interaction at high energies, but 
in many cases the RPA approximation (which amounts 
to a summation of ladder series of vertex renormalization 
diagrams in the spin channel and generally accounts for 
a Stoner-type instability at some Q) yields quite reason- 
able effective interactions, particularly near a magnetic 
instability. 

Which of the two approaches better describes FeSC is 
debatable. The value of the ordered magnetic moment 
is quite small, at least in the FESCs, and magnetic ex- 
citations measured by neutron scattering die off at ener- 
gies 100 — 150meV. This behavior is consistent with the 
idea that magnetism comes from relatively low-energy 
fermions. The one-loop parquet RG equations for the 
coupled flow of the vertices are rigorously justified only 
when the fermionic dispersion can be approximated by 
a k 2 dependence. This also does not extend too far in 
energy, i.e., the upper limit for the RG approach is also 
only a fraction of the bandwidth. 

But are these two approaches fundamentally different? 
To compare them, we remind the reader that in the RG 
treatment, the flow of magnetic and superconducting ver- 
tices remain coupled only down to energies of order Ep. 
At energies E < Ep, each vertex flows independently, 
and the flow of the pairing vertex (r sc ) becomes the 
same as in the BCS theory (aT sc /d\ogE = r 2 c . From 



this perspective, the real result of RG as far as pairing is 
concerned, is the renormalization of the pairing interac- 
tion from its bare value to the renormalized one at Ep. 
[This is indeed only true if the RG flow does not reach 
a fixed point down to Ep, but this is likely to be the 
case for FeSC simply because T c and T/v are both small 
compared to Ep ~ O.leV.] It is quite likely (although 
not guaranteed) that the angle-independent components 
of the renormalizcd interactions can be reproduced by 
choosing some other input parameters in the orbital ba- 
sis, i.e., the effect of the RG flow could be absorbed into 
the modification of the bare theory. The same is true 
for the RPA/SF approach - the angle-independent com- 
ponents of the new SF interactions likely can be repro- 
duced the renormalization of input parameters, although 
the needed values of U, V, J, J' may seem quite exotic. 

The situation is a bit more tricky for the angle de- 
pendent components of the interactions. The RG with 
coupled magnetic and superconducting vertices can be 
rigorously justified at weak coupling only if angular de- 
pendence of the vertices are (i) weak and (ii) are pre- 
served under RG. This is how the analytic RG flow has 
been obtained 12 . fRG does include some variation of the 
angle dependence of interactions, but this goes beyond 
justifiable logarithmic accuracy. Given that fRG and an- 
alytic RG yield virtually identical results, it seems as if 
the only effect of RG renormalizations down to Ep is the 
rescaling of the overall magnitudes of the interactions. 
If the effect of adding the SF component to the interac- 
tion in RPA/SF approach also predominantly gives rise 
to rescaling of the magnitudes of the interactions, then 
the two approaches are not fundamentally different. 

This does not imply that the RPA/SF and RG ap- 
proaches are equivalent. Rather, the implication is that 
the outcome of applying each of the two formalisms is the 
new "bare" theory with new input parameters. These 
new parameters do differ somewhat between RPA/SF 
and RG, but the difference should not matter much if the 
pairing symmetry and the gap structure are quite robust 
with respect to parameter variations. This is another is- 
sue that one can straightforwardly verify using the LAHA 
formalism which allows one to continuously change the 
parameters. At this stage the reader might wish to peek 
into Sec. VI (Conclusions), where the above mentioned 
questions are answered within the scope of this work. 

This paper is organized as follows. In the Sec. II we 
briefly discuss the RPA/SF and LAHA formalisms and 
outline the computational procedure. In Sec. Ill we dis- 
cuss how the gap symmetry and structure are affected 
by the SF component of the interaction. In Sec. IV 
we show the results for weakly electron-doped and hole- 
doped FeSC, and in Sec. V we discuss strongly electron- 
and hole-doped FeSC, which contain only electron or 
hole pockets, respectively. We present our conclusions 
in Sec. VI. A short summary of this work is presented in 
Ref 37 
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II. THE RPA/SF AND LAHA FORMALISMS 
AND THE COMPUTATIONAL PROCEDURE 



We first briefly describe the RPA/SF and LAHA ap- 
proaches and then outline the computational procedure. 



with 



A. RPA/SF formalism 



The approach and its application to FeSC have been 
discussed in detail in several recent publications 2 ' 5 ' 6 , so 
we will be brief. The point of departure for RPA/SF 
is a 5-orbital model with intra-orbital and inter-orbital 
hopping integrals and density-density (Hubbard) and ex- 
change intra-orbital and inter-orbital interactions, given 
by 



Hi, 
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where Cf sa is the annihilation operator for electron on 
lattice site / with orbital index s and spin er, nf s = 
n fM + n f,si, n fsa = J fscr c fsa is the density operator, 

Sf s = (l/2)cy- SQ ,(? ct( 3C/ S( 3 is the spin operator, and a = 
— a. We have designated different symbols for the intra- 
orbital Coulomb interaction U, interorbital Coulomb in- 
teraction V, interorbital exchange J and "pair hopping" 
term J' for generality, but note that if they are gen- 
erated from a single spin-rotational invariant two-body 
term they are related by J' = J/2 and V = U — 5 J/4. 
The model parameters for hopping integrals (36 total) are 
obtained from the fit to density functional theory (DFT) 
band structure 38 . 

The SF component of the interaction is now obtained 
by summing up second and higher-order ladder diagrams 
in matrix orbital formalism. The total interaction (the 
sum of the direct, first order term and SF contribution) 
is then converted from orbital to band basis by dressing 
it by matrix elements associated with the hybridization 
of five Fe orbitals. The end result of this procedure for 
the purposes of the analysis of superconductivity is the 
effective BCS-typc Hamiltonian in the band description 

U = Y e »(k)c l t k c lk + Y r y"( k > k/ )4k c l-k c ik'Ci-k' 

(2) 

The quadratic term describes low-energy excitations near 
hole and electron FSs, labeled by i and j, and the inter- 
action term describes the scattering of a pair (k t, — k I) 
on the FS i to a pair (-k' f, k' I) on the FS j. The 
effective singlet interaction Tij(k, k') is then given by 

r«(k,iO = Y <(-k)<(k)R C [r^(k,k',o)] 

s,t,p,q 

x <(k')< ho, (3) 



i[/ c x r^(k-k>)LT 



tq 



ps 



(4) 



Here the Xi PA an d Xo~ PA describe the spin-fluctuation 
contribution and orbital (charge)-fluctuation contribu- 
tion, respectively, a v v . is the matrix element connecting 
orbital p with the band Vj on FS j, and the matrices 
U c and U s contain interaction parameters from Eq. (1) 
as described in Ref. 2. We will call an approximation 
to the total interaction which includes only the two first- 
order terms in Eq. (4) as constant or non-spin-fluctuation 
(NSF), while the total interaction includes the third and 
fourth spin-fluctuation terms (SF) as well. In principle 
this interaction includes also charge/orbital fluctuations 
via U c , but as these arc negligible for realistic parame- 
ters, we use the simpler "spin fluctuation" designation. 

Throughout this paper, we consider the pairing in the 
unfolded Brillouin zone, leaving aside the issue of possible 
changes of the gap due to folding. These are particularly 
important in the case of the 122 systems with IA/mmm 



symmetry 



14.36 



In the unfolded zone, the 2D electronic 



structure of weakly and moderately electron-doped FeSC 
consists of two near-circular FSs centered at T point (k 
I i and two elliptical electron FSs centered at X and Y 
points, k = (0,7r) and k = (ir, 0), respectively. For hole- 
doped and some undoped FeSC, there exists also another, 
third, hole FS located at M point k = (ir,ir). In the 
folded zone, all three hole FSs are centered at (0, 0), and 
the two electron FSs move to (7r, it) and are hybridized 
through the coupling via a pnictogen. 

The interaction r^^k, k') contains all pairing compo- 
nents for tetragonal {D^h) lattice symmetry: A\ g (s- 
wave), Big (d x 2_ y 2), B 2s (d xy ) and A^ g (g-wave). The 
s-wave gap A s (k x , k y ) is symmetric under k x — > ±fc x , k y , 
the d x 2_ y 2 gap A ( j^ 2 _ y2 (fei k y ) changes sign under k x — )■ 
k y , and so on. We focus here on s- and d x 2_ y 2-w&ve 
symmetries. 

With these considerations in mind, the BCS gap equa- 
tion then becomes the eigenvalue problem 

r dk! i 

-Yf C] ^2^o^y r ^ k ^ k ^ A --^) - KA ^ F) 

(5) 

where a is either s or d x 2_ y 2. For a circular FS, Eq. 5 is 
simplified to 



Y fj ^A^T^VOA^WO = A Q A Q ,^) (6) 



where Np,j = mj/2ir is the density of states at the FS 
j, and ip and <p' are the angles along the FSs i and j, 
respectively. 

Eqs. 5 and 6 are integral equations which in general 
can be solved only numerically. Taking M points on each 
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FS, one obtains M cigcnfunctions and M different A's in 
each of the two pairing channels. For s-wave, some of 
cigcnfunctions correspond to an s ++ gap, while others 
correspond to an s^ gap. The eigenfunction correspond- 
ing to the largest positive X a describes the pairing state 
immediately below T c . 



LAHA formalism 



the only "memory" about d-w&ve is that the sign of a 
constant changes between the two electron FSs. 

The implication of this result is that, within LAHA, 
angle-independent and cos 20 terms must appear to- 
gether in both s-wave and d-wave components of the in- 
teractions, and with comparable magnitudes. A simple 
analysis then shows that the form of the interaction de- 
pends on whether it involves hole or electron FSs. For 
the interaction between fcrmions on a hole FS, in LAHA 



The generic idea of LAHA approximation is to model 
ry(k, k') by a rather simple function of the two mo- 
menta, such that the gap equation can be solved and ana- 
lyzed analytically. In cuprates, numerous groups approx- 
imated d x 2_ y 2 gap by the first harmonic cos k x — cosk y 
(cos 20 for a circular FS) and neglected higher harmon- 
ics like cos 60, cos 100 etc. The smallness of cos(4n + 2)0 
terms with n > 1 does not follow from any underlying 
principle, but numerically cos 20 approximation works 
rather well, at least at and above optimal doping. 

Such an approximation should generally work even bet- 
ter for FeSC because all FSs are small and even electron 
ones are almost circular. By analogy with the cuprates, 
one may try to approximate s-wave eigenfunction by a 
constant along each FS, and approximate d x 2_ y2 gaps by 
cos 20. There is a caveat, however — such approxima- 
tion is only valid for the gaps along hole FSs which are 
centered at the points along k x = ±k y [i.e., at k = (0,0) 
and (±7r, ±7r)]. Electron FSs are centered at X and Y 
points, which by itself are not k x —> ik y symmetric. As 
a result, some of s-wave gap functions, like cos k x + cos k y 
behave as ± cos 20 along electron FSs, while some of d- 
wave gap functions like cos k x — cos k y arc approximated 
by constants on the two electron FSs. In the latter case, 



r/ l / l (0,0') = A hh + A hh cos 20 cos 20' 



(7) 



where and 0' are the angles along a hole FS (measured 
relative to the k x axis), and A and A terms are s-wave 
and d-wave components, respectively. For the interaction 
between fermions from a hole and an electron FSs, 

Y eh (<t>,6) = A eh (1 + 2acos2(9)+i e/l cos20(l + 25 cos 20') 

(8) 

where 9 is the angle along an electron FS (again, mea- 
sured relative to the k x axis). Finally, for the interaction 
between fermions from an electron FS, we have in LAHA 



r ee (M') 

+A e 



A ee [1 + 2a (cos 26 + cos 26') 



4/3 cos 26» cos 26'} 
(9) 



1 + 2a (cos 26 + cos 26') + A/3 cos 26 cos 26' 

The s-wave and d-wave components look identical, but 
they transform differently between intra and inter-pockct 
interactions involving the electron FSs. 

Below we present the full LAHA result for ry(kp, k^) 
for the case when the FS consists of two hole and two 
electron pockets. The extension of the case of three hole 
FSs is straightforward. We have 



r?i 2 h 2 (0, 0') 

F/ il h 2 (0, 0') 



F/i 1 e 2 ( 
F/i 2 ei ( 
Ffa 2 e 2 ( 

Tpi pi ( 



\6') 
1,6') 
1,0') 



U hlhl + Uhihi cos 20 cos 20' 
Uh 2 h 2 + U h2 h 2 cos 20 cos 20' 
Uhti^+Uh^ cos20cos20' 

U hie (l + 2a hie cos26) + U hie (l + 2a hie cos 26) cos 20 
U hie (l - 2a hie cos26) + U hie (-1 + 2a hie cos 26) cos 20 
U h2e {l + 2a h2e cos 2(9) + U h2e (l + 2a h2e cos26») cos 20 
U h2e (l - 2a h2e cos26>) + U h2e (-1 + 2a h2e cos26) cos 20 
U ee [1 + 2a ee (cos 26 + cos 26') + A/3 ee cos 26 cos 26'} + U ee 

U ee [1 - 2a ee (cos 26 + cos 26') + A/3 ee cos 26 cos 26'} + U ee 

U ee [1 + 2a ee (cos 26 - cos 26') - 4/3 ee cos 26 cos 26'} + U ee 



(10) 



1 + 2a ee (cos 26 + cos 26') + 4(3 ee cos 26 cos 26' 
1 - 2a ee (cos26» + cos26»') + 4f3 ee cos 26 cos 26' 
-1 - 2a ee (cos26> - cos26»') + 4/3 ee cos 29 cos 29' 



Intra-pocket s-wave components U^hi and U ee rep- Coulomb (Hubbard) interaction and are positive (repul- 
rcscnt the total strength of intra- and inter-orbital sive). The signs of Uh ie are determined by the rela- 
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tive magnitudes of intra-orbital U and inter-orbital V 
terms. For a toy model of two orbitals which hybridize 
to give one hole and one electron band, Uhe > if intra- 
orbital Hubbard repulsion U exceeds inter-orbital V , and 
Uh e < for V > U 11 . It is generally expected that the 
intra-orbital interaction is the largest, and below we as- 
sume that Uhe is positive. If, however, V > U, the inter- 
action between electrons and holes is attractive in which 
case it favors a conventional s-wave superconductivity 14 . 
The signs of d-wave components of the interactions are 
determined by a more subtle balance between different 
Us and Js. 

The equation for A Qi and X a is the same as Eq. 6. 
However, now we don't have to solve integral equation 
on the gap because for Tij given by Eq. 10, s-wave and 
d-wave gaps have only four components each. For s-wave 



gap we have 

A hl 0) = A hl , A ei (9) = A e + A e cos 29 (11) 
A h2 (0) = A h2 , A e2 (9) = A e -A e cos 29 

and for d-wave gap we have 

A hl (4>) = A hl cos 20, A ei (9) = A e + A e cos 26» (12) 
A h2 (4>) = A h2 cos 20, A £2 (9) = -A e + A e cos 29 

To obtain A Q and eigcnfunctions we need to write down 
and diagonalizc 4x4 matrix gap equations. For this pur- 
pose, it is convenient to introduce bare pairing vertices 
Ar' and solve the set of conventional ladder equations 
for the full A, (0) . For the s-wave solution we obtain 



/ 1 + u hl h t L u hlh2 L 
Uh 1 h 2 L 1 + Uh 2 h 2 L 



Uh ie L 



u h 2 eL 



2uh 2 eL 
1 + 2UeeL 



2a.h 1 eUh 1 eL 
2dh 2 eUh 2 eL 

2a ee u e eL 




( A l \ 



A* 

V a°; 



(13) 



\ 2a>h 1 eUh 1 eL 2ah 2 eUh 2 eL 4:a ee U e eL 1 + 4/3 ee U ee L 

where Uij = Uij{Np,iNp,j) 1 /' 2 are dimensionless interactions, L ~ In {^f~^ C^c is measured in units of energy). The 

eigenvalue \ s and the corresponding eigenfunctions are obtained by diagonalizing this matrix equation and casting 
the result as 



(l-[AjL)(A hl , A h2 , A e , A e f 
For d-wave gap function we obtain 
(l 



(0) 



A (°>, A(°))* 



Ufi-i ho T 

u *l' T 
2 



u hih 2 T 
■ Iff h 2 h 2 



-L 



2u h2 eL 
1 + 2u PP L 



^a hl eU hl eL 
2o>h 2 eUh 2 eL 













A° 







(14) 



(15) 



Diagonalizing this set, we obtain A^ and the correspond- 
ing eigcnfunctions. 

We emphasize that the equations on X s and A^ are 4th 
order algebraic equations, hence there are effectively only 
four parameters which determine the couplings and the 
gap structure in each channel. This is the same number 
as the minimal number of interaction parameters in the 
original orbital model. 



If the two hole FSs can be treated as equal (i.e., Uhihj = 
Uhhi u h ie = u h2 e = u he and a hie = ot h2 e = a, which 
actually is quite consistent with the fits to RPA/SF, see 
below), then A^ = Ah 2 = Ah and the gap equation 
reduces to 3 x 3 set which can be very easily analyzed 
analytically. The equation for the s-wave gap in this 
case is 



1 + 2u h hL 2u he L 2a he UheL 

2u he L l + 2u ee L 2(XeeUeeU I I A e I = I A^ I (16) 



ictheUheL 4a 

ee^ee 

L 1 + 4/3 ee it ee L 
For the d-wave gap, we have 

1 + U h hL 2u he L 2ckheUUheL 



A h \ 


K 


A e 


= A° 


A e J 





A h \ 




A e 


■(3 


A e y 





U he L \ + 2UeeL 2a ee UeeL ) | A e | = | A° ] (17) 

2& he U he L ActeeUeeL 1+4/3 
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For the case of three hole FSs, we have to introduce 
three different A^ . The gap equations in s-wave and d- 
wavc channels become 5x5 sets. Still, they can be very 
easily analyzed. 

C. General Considerations 

Before we proceed further, it is instructive to take a 
more careful look at the 3x3 sets Eqs. 16 and 17 to illus- 
trate issues which we outlined in the introduction. Con- 
sider first the s-wave gap equation. Suppose momentarily 
that all interactions are angle- independent, i.e., a/ le , a ee 
and p ee = 0. Then we have three solutions 

[0 

A s = ~ \ Z( M f a ) 2 +"L (18) 

{ -^^ + ^^) 2 +ul e 

A positive (attractive) A s emerges only when u^ h > 
u ee Uhhi i.e., when inter-pocket pair hopping interaction 
term exceeds intra-pocket repulsion. The eigenfunction 
corresponding to the positive X s is a sign-changing 
gap: A e = — A/,,, the one corresponding to the negative 
A s is a conventional s ++ gap: A e = A/ t . 

Suppose next that the angular dependence of the in- 
teraction is present. Then ± cos 29 components of the 
gaps along electron FSs become non-zero. Solving the 
cubic equation for X s , we now find three non-zero solu- 
tions. When u 2 eh > u ee Uhhi the solution with the largest 
A s gradually evolves from the one which already existed 
for constant interactions. When ahe, ot eei and /3 ee in- 
crease, the cos 29 component of A e grows, and at some 
point gets larger than A e , and the gap develop nodes. 
This is one scenario. Another one comes from the analy- 
sis of the region u 2 eh < u ee Uhh, where no solution was 
possible without angular dependence of the interactions. 
In that case, one of A s was zero. When the angular de- 
pendence is included, this X s becomes non-zero, and its 
sign is determined by the sign of 

S = UeeUhh {ale ~ Pee) + u\ h (p? eh + p ee ~ 2a ee a eh ) 

(19) 

When S > 0, A s > 0, i.e., the system develops an attrac- 
tion in s-wave channel exclusively because of the angular 
dependence of the interactions. The eigenfunction cor- 
responding to such A s has nodes even when the angular 
dependence of the interactions are weak. The physics pic- 
ture is that the system finds a way to minimize the effect 
of strong intra-pocket repulsion u ee by inflating cos 29 
components of the gaps along the electron FSs, because 
these components do not couple to angle-independent 
component of the interaction. 

Note that the sign of S is predominantly determined 
by the interplay between the angular dependence of 
electron-hole and electron-electron interactions. When 
a e h 3> a ee , or when a\ e > p ee , S > even when 
UeeUhh 3> igj,, i.e., superconductivity with nodes 



develops despite the fact that intra-pocket repulsion is 
the strongest. In particular, S is definitely positive 
if only electron-hole interaction has momentum depen- 
dence, i.e., if a ee = Pee = 0. In this situation, S = 
u lh a lh > o. 

Consider next the limit when hole FSs are absent and 
only electron ones are present. At the first thought, s- 
wave pairing is impossible. On a more careful look, how- 
ever, we find from Eq. (16) that, even if we set u e h = 0, 
one of A s for s-wave pairing is still positive if p ee < Qg e , 
no matter how small angular dependence of electron- 
electron interaction is. The eigenfunction for this solu- 
tion again has nodes. The physics reasoning is the same 
as in the case we just considered: the angle-independent 
part of the interaction is repulsive, but the system finds a 
way to overcome this strong repulsion by inflating cos 29 
components of the gaps along the electron FSs. 

We next turn to the d-wave gap equation ( Eq. 17). 
The generic reasoning parallels the one for s-wave case, 
namely for interactions independent on cos 29 attractive 
d-wave solution exists when u\ h > u ee Uhh, and for angle- 
dependent interaction one of Xd is positive (attractive) 
even if this condition is not satisfied, but S, which is a 
d-w&ve analog of S from Eq. (19), is positive. There is, 
however, one crucial distinction with the s-wave case: the 
d-wave interactions Uhh an d u ee are not necessary posi- 
tive. In particular, u ee is the difference between angle- 
independent components of intra-pocket and inter-pocket 
interactions between electron pockets. Once the inter- 
pocket interaction is larger (e.g., when magnetic fluctu- 
ations are peaked at (tt, vt)), u ee < 0, and the system 
develops an attraction in the d-wave channel, even when 
d-wave electron-hole interaction is weak. This is partic- 
ularly relevant for the case when only electron FSs are 
present, i.e., within BCS approximation u e h can be set 
to zero. If u ee < in this case, d-wave solution emerges, 
with the sign-changing gap on the two electron FSs. An 
alternative possibility is that u ee > 0, but a\ e > /3 ee , and 
d-wave attraction is produced by angle-dependent part of 
d-wavc interaction. In this situation, cos 29 component 
of the d-wave gap is large, and the gap has nodes on the 
two electron FSs. 



D. Computational Procedure 

We use Yij ;(kp, k^) obtained in the RPA/SF approach 
as inputs and fit their functional forms by Eq. 10. This 
gives us Uij, c*ij and Pij. We assume for simplicity that 
Np t j are the same for all FSs, convert Uij into Uij = 
NpUij and analyze 4x4 and 5x5 gap equations. We 
then vary parameters and check how robust the solutions 
arc. The approach can be extended to the case when Npj 
depends on j, but this dependence very likely does not 
change the physics. 

We analyzed six different sets of parameters and sev- 
eral different doping levels which correspond to either 
electron-doping or hole doping. The parameters are pre- 
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scntcd in Table I. The results for all sets of parameters 
are quite similar, and below we show the results only for 
representative cases. We also analyzed the case of large 
electron doping, when there are no hole FSs. 

TABLE I: Six sets of parameters used in comparison of 
RPA/SF and LAHA. We used different values of the chemi- 
cal potential /i for each set (fj, = —0.30 to +0.30), the range 
covers hole doping(negative) to electron doping (positive). 



set # U J V_ 

set 1 1.67 0.21 1.46 

set 2 1.0 0.25 0.69 

set 3 1.2 0.0 1.2 

set 4 1.0 0.9 -1.25 

set 5 1.0 4.0 -4.0 

set 6 1.0 0.9 -0.7 



The use of LAHA is meaningful only if the fit of the ac- 
tual TijCkp, k' F ) by Eq. 10 is accurate enough. In Figs. 1 
and 2 we show fits for the representative case of /i = 0.08 
(electron doping with n e = 6.12, while undoped case is 
n e = 6) and interactions from the set #1 (U — 1.67, 
J = J' = 0.21, V = 1.41) considered in Ref. 1. We 
set one of momenta to be cither along x or along y axis 
on one of the four FSs and vary the other momentum 
along all four FSs. The parameters for s-wave and d- 
wave interactions extracted from the fit are presented in 
Table II. 

The original (kp , k' F ) are in the units of eV. Within 
our approximation of an equal Nfj on all FSs, the di- 
mensionless Uij are r^kf, k F ) times the rescaling fac- 
tor. This factor affects the overall scale of the eigenvalues 
A, but doesn't affect the sign and relative magnitudes of 
A s and A<2, which we need. For simplicity, we then set 
this rescaling factor to be equal to one. 

We see that the fits are quite good. For the case of no 
SF term (NSF), all T tJ (k F , k' F ) are fitted well by Eq. 10. 
For the full interaction (i.e. SF included), all (kp, k F ) 
are fitted well, the small discrepancies can be cured by 
including cos4# terms into the interactions. 

The fits for other parameters arc quite similar. For 
example in Figs. 7-9 below we show the "best" and 
the "worst" fits of interactions for the set #2 for three 
different values of \i. For positive fi, there are 4 FSs, and 
for negative /i, there is an additional, 5th hole FS. The 
fits are not perfect, but are quite good for all practical 
purposes. 

We continue below with the set #1 for the discussion 
on how SF contribution affects the pairing symmetry and 
the gap structure. The trend is the same for all other 
sets. And just for a change we will look at set #2 for the 
discussion of how gap symmetry and structure change 
with doping. We do so to limit the number of figures 
and not overwhelm the readers. We will point out in the 
text if differences arise. 



III. SENSITIVITY OF THE GAP STRUCTURE 
TO THE FORM OF r/yfkfjk'j.) 

We first use the parameters extracted from the fit for 
the set #1 and solve the 4x4 gap equation within LAHA. 
The results for the case of no SF component are shown 
in Fig. 3. For comparison, the gaps obtained by the full 
numerical solution within RPA are also presented. The 
LAHA and RPA solutions are almost equivalent, which 
is another indication that the LAHA fit works quite well. 
The results for the full interaction, with SF component, 
are also presented in Fig. 3 along with the numerical so- 
lution for s-wave gap within RPA/SF. The LAHA and 
RPA solutions again agree very well. The only real dif- 
ference between RPA/SF and LAHA is that the value 
of the gap along the outer hole FS is somewhat larger 
in RPA/SF. Like we said, this is the consequence of the 
fact that to fit r^ 2ie (</), 9) in LAHA one needs to add 
cosAB components. Once we include these components, 
the gap along the outer hole FS goes up, bringing the 
LAHA result even closer to RPA/SF. 

We see therefore that LAHA approximation works 
quite well both for the bare interaction and the full inter- 
action with SF contribution. We verified that the near 
equivalence between gaps obtained within RPA/SF and 
LAHA holds for all other sets of parameters from Ta- 
ble I. This gives confidence that the physics can be 
understood by analyzing 4x4 s-wave and d-wave gap 
equations within LAHA. We begin with s-wave case. 

A. s-wave case 

1. Modification of Tij(]iF ,k' F ) by spin fluctuations 

We first take a closer look at Table II. Comparing the 
values of u hlhl ,u h2h2 ,u hl h2,u hl e,Uh 2 e, and u ee we see 
two trends. First, once SF contribution is added, there is 
overall enhancement, roughly by a factor of three, for all 
interactions, including interactions within hole pockets. 
On top of this overall enhancement, there is another ef- 
fect - electron-hole interactions Uhe further increase com- 
pared to Uhh and u ee . This additional increase is by fac- 
tor of 3 — 4, such that the total increase of Uh ie is by a 
factor of 10 [the increase of Uh 2 e is a bit smaller]. 

We attribute the overall increase to the "bare" SF 
interaction term u 2 xo(q) and the additional increase 
to the relative enhancement of the RPA susceptibility 
x(Q)/xo(Q) near the momentum transfer Q — (it, 0) or 
(0, 7r). The reasoning is based on the comparison with the 
results of Ref. 2 for the bare and RPA-renormalized spin 
susceptibility. First, the product UijXoO^-i ~ kj) remains 
roughly constant in k— space, if we use «y from Table II. 
Second, the RPA-renormalized x(Q) is 3-4 times larger 
than Xo(Q), while x(0) « Xo(fJ). (Ref. 2). This logic also 
applies to the interaction between electron pockets, for 
which u ee is the average between intra-pocket interaction, 
for which \ = x(0) ~ Xo(0), and inter-pocket interaction 
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FIG. 1: Fits of the actual interactions Fij(\iF,\c' p ) by LAHA. Symbols represent interactions computed numerically for the 
5-orbital model using LDA band structure, solid lines are fits using Eq. 10. The fit is for the set #1 for bare interactions (no 
SF component). The chemical potential is /j, — 0.08 which corresponds to electron doping. We set in Fij(k.F,]i' F ) to be 
either along x or along y direction on a given FS (its location is specified in the title on top of each figure) and varied k^- along 
each of FSs. The angle 6' is measured relative to k x . The fit is amazingly good. 




Angle/(2n) Angle/(2n) Angle/(2n) Angle/(27i) 



FIG. 2: The same as in Fig. 1 but for the full interaction, with SF component. The fit is again amazingly good. We verified 
that the slight discrepancies are removed if in LAHA we add a 40 harmonics to the interaction. 



for which \ = x( 7r : 7r ) ~ 2.4xo(7r, t)- The total increase 
of u ee is then expected to be 1 + 2((1 + 2.4)/2) = 4.4, 
and we see from the Table II that u ee increases by quite 
similar factor of 4. 

Compare next the angular parts a^e, ®ee and f3 ee . We 
see from Table II that cthe and a ee do not change much. 
The term /3 ee does change and becomes 2.5 times smaller 
in the presence of the SF component. However, we will 
show in the next section that the gap structure is insensi- 
tive to the change of /3 ee and does not change much even 
if we set (3 ee = (see Fig. 5). 

Neglecting the change of /3 ee , we conclude that the 
SF contribution to the pairing interaction increases the 
overall magnitude of Tij and additionally increases the 
magnitudes of electron-hole interactions (uh ie and Uh 2 e 
terms) and, to lesser extent, of electron-electron interac- 
tion (the u ee term), but doesn't substantially modify the 
relevant angular dependence of the electron-hole interac- 
tion. The overall increase of the pairing interaction does 
not affect the gap structure, hence the only true effect of 
the SF term is the increase of electron-hole interactions 
compared to hole-hole and electron-electron interactions. 

As we said, controlled RG flow of the couplings gives 
rise to exactly the same effect - angular dependence are 
preserved during the flow, but the relative magnitude of 
electron-hole interaction increases. From this perspec- 
tive, RG and SF approaches, although formally different, 
describe very similar physics. 



2. Effect of the angular dependencies of electron-hole and 
electron- electron interactions 

We now analyze explicitly how sensitive is the gap 
structure to various angular-dependent components of 
Tij. Within LAHA, we can easily change angular de- 
pendence of the interactions and check the consequences. 
First, we verify how sensitive is the solution for the gap 
to the change of the angular component of electron-hole 
interaction. 

In Figs. 4 and 5 we show the results for the gap for 
different values of ah ie - Comparing these figures with 
Fig. 3, we see that the effect of o^e on the gap is dif- 
ferent for bare and full interaction. For the bare interac- 
tion, originally there is no s-wave solution with positive 
A s , but it appears once we increase a^e above a cer- 
tain threshold. This can be easily understood by analyz- 
ing 4x4 gap equation: for bare interaction intra-pocket 
repulsions u^h- and u ee are stronger than inter-pocket 
Uh ie , hence the pairing can only be induced by angle- 
dependent components of the interaction, when the fac- 
tor S, given by Eq. (19), becomes positive. For original 
parameters S < 0, but once we increase a>h ie , S eventu- 
ally changes sign and the solution with X s > appears. 
The gap function for this induced solution has strong os- 
cillations along electron FSs and has accidental nodes. 

Consider next the full interaction with SF component. 
Uh ie are now enhanced, and the solution with X s > 
exists even if we set a,h ie = (see Fig. 5). The only dif- 
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TABLE II: s- and d- wave parameters for the set #1 with /i = 0.08. Here and henceforth NSF and SF mean the bare interaction 
without SF component and the full interaction with SF component, respectively. 



s— wave 




u h 2 h 2 








Uh 2 e 


0/12C 




Otee 


Pee 


NSF 


0.8 


0.76 


0.78 


0.46 


-0.24 


0.4 


-0.30 


0.77 


0.14 


0.09 


SF 


2.27 


2.13 


2.22 


4.65 


-0.34 


2.29 


-0.22 


3.67 


0.15 


0.04 


d-wave 










&h ie 








&ee 


Pee 


NSF 


0.7 


0.66 


-0.68 


-0.25 


-0.58 


0.24 


-0.42 


0.11 


-0.5 


0.25 


SF 


1.50 


1.40 


-1.50 


-3.73 


-0.44 


1.44 


-0.32 


1.03 


-0.49 


-0.02 




FIG. 3: . Top: s- and d— wave solutions obtained by applying LAHA for the gap structure (left pair) for the set #1 (NSF) with 
(j, — 0.08. s- and d— wave solutions obtained numerically from the 5-orbital model based on the LDA band structure(right pair) 
are also shown. There are no solutions with positive A, so we show the solutions with the smallest negative A in both channels. 
The agreement of the gap structure with LAHA results is quite good. For LAHA, the s-wave solution has A s = —1.02 and 
d-wave solution has Ad = —0.99 (A s /Ad ~ 1.0). A s and Ad obtained from the numerical approach are negative, and their ratio 
is A s /A,j = 0.6. Bottom: The same, but now with the SF component of interaction included. The couplings are now positive: 
A a = 4.6 and X d = 4.8 (A s /A d « 0.96) for LAHA and » 1.1 for RPA/SF calculations for the 5-orbital model 1 . 



ference between the solutions with small and larger ah ie 
is that in the first case gap has no nodes on the two 
electron FSs. We see therefore that for full interaction 
the role of oi] lie is merely to modify already existing 
solution and add angular variation to the gap along the 
two electron FSs. 

We next analyze how sensitive is the solution for the 
gap to the change of the angular component of electron- 
electron interaction (a ee and /3 ee terms). In Fig. 5, bot- 
tom panel, we show the results for the gaps obtained with 
a ee and f3 ee increased by two, reduced by two, and set to 
zero. Comparing this figure with Fig. 3(bottom panel) 
we see that the changes in a ee and f3 ee lead to very lit- 
tle changes in the gap structure (and the As). The gap 
remains very much the same as in Fig. 3, even if we set 
a ee and /3 ee to zero. The implication is that, for full in- 
teraction, the gap structure is determined by the angular 



dependence of electron-hole interaction, while electron- 
electron interactions can be we approximated by angle- 
independent u ee term. 



Nodal vs non-nodal s-wave gap for 4 FSs 

Previous studies of the gap structure for 4 FSs within 
RPA/SF formalism yielded the gap with accidental nodes 
along the electron FS 1,2 . The solution without nodes 
only appears for an electronic structure which contains 
an additional, 5th hole FS centered at (ir, tt). The re- 
sults of the previous subsection imply that this result is 
the likely outcome, but is non-universal: for the full in- 
teraction the gap can be either nodclcss or with nodes 
depending on the magnitude of a^e and on the interplay 
between u\ e and u ee Uhh (for simplicity, we set two hole 
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FIG. 4: Effect of electron-hole interactions: s-wave solutions for the set #1 (NSF) for ft = 0.08 with the angular parts of 
electron-hole interactions {pth^e) set to zero and increased by a factor of 2, respectively (left and center). The values of \ a are 
-0.16 and 0.226, positive \ s is for enhanced a^e. In the right figure, we introduce additional enhancement of the overall factors 
of the electron-hole interactions by ~ 4 and this causes nodes in the gap to disappear. 
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FIG. 5: The dependence of the gap structure for s-wave solution on the strength of angle-dependent parts of the interactions. 
We used set #1(SF) for fi — 0.08. The solution for the original ceh ie is shown in Fig. 3. Top panel: the gap structure for the 
cases when the angular parts of electron-hole interaction (a^e) are increased by factor of 2, reduced by factor of 2, and set 
to zero, while the angular parts of electron-electron interaction (a ee and /3 ee ) are kept at original values. The values of X s are 
7.3, 3.4, and 2.5 respectively. The angular dependence do change, and, in particular, the nodes along the hole FS disappear 
when we set ah ie = 0. Bottom panel: the same but now we keep aj, ie at original value and change the angular parts of the 
electron-electron interactions (a ee and /3 ee terms): increase them by 2, reduce by 2, and set to zero, The A's are 4.5,4.7 and 4.0 
respectively. We see that this has very little effect on the structure of the gap. 



FSs to be equal). When uf le < (u ee u hh ), the only possi- 
bility is s^ 1 pairing with nodes along electron FSs, in the 
opposite case the solution can be either nodal or nodeless 
depending on ate- The larger u\ e / (u ee Uhh) is, the larger 
ate one needs to make SC nodeless. 

To further illustrate this, in Fig. 4(right) we show the 
solution of 4 x 4 set for the gap for the case when we 
increase Uh ie and Uh 2e keeping all other parameters, in- 
cluding a^e, hxed. We see that the nodes on electron 
FSs disappear for substantially large Uh ie and Uh 2 e- 

In RPA/SF analysis, a way to further increase ui lie 
compared to other parameters is to bring the system 
even closer to AFM instability. For the parameter set 



we are dealing with, one needs to increase Uh ie by quite 
substantial amount, so a fine tuning to AFM instability 
is required to eliminate the nodes. But still, no-nodal 
solution definitely exists in some range of parameters. 



B. d-wave case 

The consideration for the d-wave case proceeds simi- 
lar to the s-wave case. The d-wave pairing again can be 
due to attraction coming from angle-independent parts of 
the interaction, or it can be induced by angle-dependent 
parts of the interactions, even when < in the absence 
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of anglc-dcpcndent terms. We refrain from discussing all 
cases as in many respects the analysis for e?-wave pairing 
parallels the one for the s-wave case. We point out, how- 
ever, that the negative Xd for the case of bare interaction 
(NSF) is much smaller in magnitude than the coupling in 
the s-wave channel, i.e., there is "less repulsion" in the d- 
wave channel. As a result, it takes less to convert d-wavc 
repulsion into attraction. This can be achieved by shift- 
ing the values &h ie , fire, and /3 ee . Once SF piece is added 
to the interaction all <i-wave components increase, but, 
just as for s-wave case, the largest increase is for electron- 
hole interaction. This increase is large enough such that 
the gap equation develops a solution with > even if 
we set angular components of electron-hole and electron- 
electron interactions to zero. In the latter case, the gaps 
along the two electron FSs are ±A e . When we don't set 
angle-dependent components to zero, electron gaps ac- 
quire cos 29 components, but, as we see in Fig. 3, this 
component is quite small. 

We also note that the values of A s and Xd for the full 
interaction are quite close. For the set #1 which we are 
discussing in this section, we found that, within LAHA, 
Ad is actually a bit larger than X s , for the parameters 
extracted from the fit, but A^ and A s are very close and 
A s becomes larger already after a small change of param- 
eters. The solution for the gap within RPA also shows 
that, for this set, A s and Xd are very close, but it yields 
A s > Xd- In any event, however, s-wave coupling A s def- 
initely becomes stronger than Xd once the system comes 
close enough to an antiferromagnctic instability. 



IV. DOPING EVOLUTION OF s-WAVE AND 
d-WAVE GAP FUNCTIONS 



To understand the evolution with doping, we choose 
set #2 for definitcness and consider the gap structure for 
positive and negative values of the chemical potential. 
Positive /i corresponds to electron doping and negative 
jj, correspond to hole doping. For positive fi, electronic 
structure consists of two hole and two electron pockets, 
and the size of hole pockets gets smaller as /x increases. 
For negative /i, the electronic structure contains an ad- 
ditional, 5th hole pocket, centered at (tt,tv). As before 
we shall denote the electron pocket at (tt,0) as e%, the 
electron pocket at (0, 7r) as e2, the inner hole pocket at 
(0, 0) as h\, the outer hole pocket at (0, 0) as h 2 and the 
hole pocket at (tt,tt) when present as /13. In this section 
we only consider moderate doping, when both hole and 
electron FSs are present. We consider the limiting case 
of large electron and hole dopings in the next section. 
For simplicity, we only present results for full interaction 
with the SF component. 



A. Electron doping 

We considered several values of electron doping (posi- 
tive n). Below we present the results for the representa- 
tive case of n = 0.05 (n e = 6.09). The FS is presented 
in Fig. 6 In Table III we present LAHA s-wave and d- 
wave interaction components from Eq. (10), obtained by 
fitting RPA/SF results. The actual fits are presented in 
Fig. 7, where we present the "best" and the "worst" fits. 
We see that even the "worst" fits are actually quite rea- 
sonable. In the same figure we show s-wave and d-wavc 
gaps corresponding to the largest A s and A^. 

The tables show the same trends that we discussed in 
the previous section, the most relevant one is the overall 
increase of all interactions when SF component is added, 
and further increase of electron-hole interaction in both 
s-wave and d-wave channels. In this Section we focus on 
the features associated with the dependence on doping. 
For briefness, we show only the case of full interaction 
with SF component. 

First, we see from Fig. 7 that for positive fi (when 
the electronic structure consists of 4FSs) , s-wave solution 
with As > has nodes on the electron FSs. The d-wave 
solution with A^ > has symmetry-related nodes on the 
hole FSs, but no nodes on electron FSs. That > is 
the combination of the two effects: (i) the d-wave compo- 
nent of the interaction between electron pockets is neg- 
ative (u ee < 0), i.e., there is a direct d-wave attraction 
between electron pockets, and (ii) d-wave components of 
electron-hole interaction Uh ie and Uh 2 e are quite large. 
Both effects give rise to > even if we set all cos 29 
components of the interactions to zero. Within LAHA, 
we can check the relative importance of the two effects by 
artificially setting one of them to zero. We see from Ta- 
bles III and IV that u ee is attractive, but very weak. The 
primary reason for A^ > are large values of electron- 
hole interaction Uh ie and Uh 2 e- Accordingly, the driving 
force for d-wave attraction is strong d-wave component 
of the pair-hopping between electron and hole pockets. 
In this respect, the mechanism is quite similar to that 
for sign-changing s-wave gap. 

How far in electron doping does this mechanism remain 
the leading one? To analyze this, we also considered the 
case of a larger /j, = 0.18 (n e = 6.23), when hole pock- 
ets almost disappear. The FS for this case is shown in 
Fig. 6, the parameters extracted from the fit are shown 
in Table IV, and the fits and the gaps are presented in 
Fig. 8. A somewhat surprising result is that there is very 
little change compared to the case of smaller /j, = 0.05, 
when hole pockets are much larger. Still, the dominant 
interaction in the (i-wave channel is the pair hopping be- 
tween hole and electron FSs. The d-wave component of 
electron-electron interaction is negative (i.e., attractive), 
but it remains very small. 

Note also that for both yu's Xd > X s , and the differ- 
ence increases as electron doping increases. This does 
not necessarily mean that d-wave is the leading insta- 
bility because T c for s-wave and d-wave superconducting 
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FIG. 6: (Top: left to right) FSs for increasing electron doping with /j, = 0.05, 0.18, 0.30 respectively, /j, — 0.18 represents the 
electron doping at which the hole FS almost disappear, while = 0.30 represents extreme electron doping where the hole FS 
completely disappear. (Bottom-left to right) FSs for increasing hole doping with /i = —0.05, — 0.18, — 0.30 respectively. Here 
electron FSs almost disappear at /j, — —0.18 and completely disappear at /i = —0.30 



TABLE III: s- and d-wave parameters for the set #2 with /j, = 0.05 



s-wave 




Uh 2 h 2 


Uh 1 h 2 




Cth t e 


Uh 2 e 


Oth 2 e 


Uee 


Ctee 


Pee 


NSF 


0.51 


0.50 


0.50 


0.32 


-0.19 


0.28 


-0.24 


0.50 


0.11 


0.08 


SF 


0.80 


0.79 


0.79 


0.79 


-0.19 


0.67 


-0.19 


0.91 


0.05 


0.05 


d-wave 


Uh 1 h 1 


Uh 2 h 2 


Uh x h 2 


Uh x e 




Ufi 2 e 


&h 2 e 


Uee 


Oiee 


Pee 


NSF 


0.38 


0.36 


-0.37 


-0.14 


-0.60 


0.14 


-0.60 


0.05 


-0.6 


0.35 


SF 


0.50 


0.49 


-0.50 


-0.39 


-0.46 


0.30 


-0.47 


-0.04 


1.5 


-0.69 



TABLE IV: s- and d-wave parameters for the set #2 with /i = 0.18. 



s-wave 




Ufi 2 h 2 


u hj_h 2 


Uh ± e 




Uh 2 e 


OLh 2 e 


Uee 


Ctee 


Pee 


NSF 


0.56 


0.55 


0.56 


0.30 


-.23 


0.29 


-0.23 


0.52 


0.12 


0.06 


SF 


0.75 


0.75 


0.74 


0.67 


-0.19 


0.67 


-0.20 


0.88 


0.10 


0.05 


d-wave 


Uhxh-i 


Uh 2 h 2 


Uhih 2 


Uh x e 


&hxe 


Uh 2 e 


&h 2 e 


^ee 


&ee 


Pee 


NSF 


0.43 


0.43 


-0.43 


-0.14 


-0.62 


0.14 


-0.62 


0.04 


-0.625 


0.44 


SF 


0.51 


0.51 


-0.51 


-0.32 


-0.50 


0.32 


-0.50 


-0.05 


0.9 


-0.6 



TABLE V: s- and d-wave parameters for the set #2 with /j, — —0.05. 



s-wave 


Uh x hx 


Uh 2 h 2 


Uh 3 h 3 


Uh ± h 2 


U hl h 3 


u h 2 h^ 


Uh ± e 


Oth t e 


Uh 2 e 




Uh^e 


Oth 3 e 


Uee 


Ctee 


Pee 


NSF 


0.50 


0.49 


1.00 


0.49 


0.20 


0.13 


0.34 


-0.18 


0.30 


-0.25 


0.61 


0.26 


0.49 


0.11 


0.10 


SF 


0.86 


0.96 


1.83 


0.89 


0.45 


0.32 


0.92 


-0.18 


0.79 


-0.21 


1.5 


0.21 


1.00 


0.11 


0.08 


d-wave 




u h 2 h 2 


Uh 3 h 3 


Uh ± h 2 


U hl h 3 


u h 2 h 3 


Uh ± e 


&h ie 


Uh 2 e 


&h 2 e 


Uh^e 


&h 3 e 


Uee 


&ee 


Pee 


NSF 


0.36 


0.36 





-0.36 





-0 


-0.15 


-0.58 


0.15 


-0.58 


-0 


-0.58 


0.06 


-0.58 


0.33 


SF 


0.51 


0.61 


0.01 


-0.56 


0.00 


-0.01 


-0.45 


-0.48 


0.39 


-0.43 


0.02 


0.90 


0.07 


-1.00 


0.46 
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TABLE VI: s- and d-wave parameters for /i = —0.18. For technical reasons we used U = 0.9, J = 0, and V = 0.9. 



s-wave 






















^fr.3 e 










NSF 


0.37 


0.37 


0.74 


0.36 


0.04 


0.10 


0.40 


0.0 


0.40 


0.0 


0.04 


0.0 


0.44 


0.0 


0.0 


SF 


0.75 


2.02 


17.2 


0.98 


-0.08 


0.41 


1.36 


0.08 


2.86 


0.02 


0.31 


-0.01 


1.40 


0.01 


0.06 


d-wave 






Uh 3 h 3 


Uhxh.2 


U hl h 3 


Uh 2 h 3 








&h 2 e 


^fr.3 e 








/3 ee 


NSF 


0.36 


0.36 





-0.36 


-0.04 


0.04 


-0.40 


-0.0 


0.40 


-0.0 


0.04 


-0.0 


0.44 


-0.0 


0.0 


SF 


0.70 


1.94 


13.6 


-0.94 


0.04 


0.33 


-1.32 


0.0 


2.85 


0.02 


0.26 


0.02 


1.45 


0.01 


0.04 




FIG. 7: Top panel. The 'best' and 'worst' LAHA fits of the interactions F(k.F,k^) obtained from RPA/SF calculations for 
parameter set #2(SF) for /j, = 0.05. Bottom panel: Gap-functions in s- and d-channels obtained in the LAHA. X s = 0.25, 
Ad = 0.37. Note that Xd is somewhat larger. 



instabilities have different prefactors. Still, a larger value 
of Xd implies that <i-wave superconductivity is certainly 
a possibility in electron-doped pnictides. A more exotic 
mixed s + id state is also quite possible, 39 but to study it 
one obviously needs to solve a non-linear gap equation, 
which is beyond the scope of this work. 



B. Hole doping 

We next consider representative cases of hole doping 
by setting fi to negative values, /j = —0.05 (n e = 5.95) 
and (i = -0.18 (n e = 5.53). In Fig. 6 we show the FS 
for these two fj,'s. The FS now has an additional hole 
pocket centered at (tt,tt). For /i = —0.05 hole and elec- 
tron pockets are of comparable size, for \i = —0.18, the 
electron pockets almost disappear. The fits to LAHA for 
/i = —0.05 are presented in Fig. 9 and the parameters ex- 
tracted from the fits are summarized in Table V. We see 
from Fig. 9 that for hole doping the situation is different 



in two aspects. First, the s-wave gap has no nodes; sec- 
ond, A s is substantially larger than Xd, i.e., s-wave super- 
conductivity is the most likely scenario. To understand 
these differences, compare Table III with Table V). We 
see that the interactions between two hole FSs at (0, 0) 
and two electron FSs do not change substantially between 
electron-doped and hole-doped cases, but for hole-doped 
case there appear additional hole-hole and hole-electron 
interactions associated with the fifth hole pocket. These 
additional interactions are weak and irrelevant for the 
d-wave component of r,;j(kp, k' F ), but are quite strong 
for the s-wave component. We analyzed the 5x5 gap 
equation for the s-wave gap, and found that these addi- 
tional interactions effectively increase angle-independent 
component of electron-hole interaction, which favors no- 
nodal s gap. As the consequence, the system develops 
ans ± solution with relatively large X s > 0, even if all in- 
teractions are set to be angle independent. When we in- 
clude the angle-dependent components, the gap acquires 
some momentum dependence along electron FSs, but for 
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FIG. 8: The same as in Fig. 7, but for different fi — 0.18. A s = 0.21, Xd = 0.35. Again, A^ is larger, and the difference between 
Xd and A s is larger than for /j, — 0.05. 




FIG. 9: The same as in Fig. 7, but for negative /i = —0.05, when extra hole FS appears. A s = 0.58, Ad = 0.31. Observe that 
now A 3 > Xd and s gap has no nodes. 



given parameters this dependence is weak and the gap 
remains nodeless. There is no "theorem" , however, that 
the s ± solution is always nodeless for five FS pockets. In 
Fig. 11 we show the gap for the full interaction for the 
set #3. We see that the gap does in fact have nodes on 
the electron FSs. 

We also note that, for the c?-wave case, the electron- 



electron interaction u ee is repulsive, and (i-wave solution 
with Xd > is again the result of relatively strong d-wave 
electron-hole interactions Uh ie & n d Uh 2e involving the two 
hole pockets centered at (0,0). 

The outcome of this study is quite simple - no- nodal 
solution is the leading instability if the angle- independent 
part of the electron- hole interaction Uhe (either the direct 
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FIG. 10: The same as in Fig. 9, but for jj, 
still, As > Ad and s gap has no nodes. 



-0.18 and U = 0.9, J = 0.0, V = 0.9 We obtain A s = 1.8, A d = 1.2. Observe that, 
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FIG. 11: s gap structure obtained within LAHA for the 
parameter set #3(SF) for fi = —0.05. The coupling A 3 = 0.79. 
The gap has nodes on electron FS despite that there is the 
third Fermi surface at (ji, ir). This nodal solution has a strong 
d-wave competitor for which A^ = 0.8. 



one or the effective one, in case of 5 FSs) is sufficiently 
large compared to (uhhUee) 1 ^ 2 ■ In the opposite case, the 
s-wave solution has nodes on electron FSs, and d x 2_ y 2 



pairing is a strong competitor. 

Consider next how far in doping the nodeless s ± gap 
remains the leading instability. To analyze this, we turn 
to the case of /i = —0.18, when the electron pockets 
almost disappear. The FS for this case is shown in Fig. 6, 
the parameters extracted from the fit are shown in Table 
VI, and the fits and the gaps are presented in Fig. 10. 

On analyzing the structure of the interactions and the 
gaps and comparing them to the case of = —0.05 we see 
that the key features survive despite the small size of the 
electron pockets. Namely, the leading instability remains 
s ± , the gap has no nodes, and the driving force for the 
pairing is the interaction between hole and (still existing) 
electron pockets. The d-wave channel is a competitor 
(\d > 0), but still, s ± state has larger A. 

The conclusion of this subsection is that, as long as 
both hole and electron pockets are present, the pairing 
instability is essentially driven by electron-hole interac- 
tion. This should obviously change at even larger hole or 
electron dopings, when only one type of pockets remain 
and the pairing (if it exists) should come from the inter- 
action either between hole pockets or between electron 
pockets. 



V. OVERDOPING 

Finally, we consider the case of strong electron or hole 
doping, when hole FSs disappear and only electron ones 
at (0, 7r) and (ir, 0) remain, and the case of strong hole 
doping, when electron FSs disappear and only hole FSs 
at (0, 0) and (tt, tt) remain. We consider the two limits 
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separately. 



A. Strong electron doping 

Strongly electron doped FeSC include recently discov- 
ered FeSe superconductors AFe2Se2 (A = K, Rb, Cs) 32 . 
T c in these materials is quite high and reaches almost 
40A'. The electronic structure of these materials in the 
folded Brillouin zone is a bit more involved because (i) 
AFe 2 Se2 has body-centered tetragonal structure what 
makes the folding of the electron FSs a more complex 
procedure than just the mixing of the two ellipses which 
would be the case for a simple tetragonal structure 14 ' 36 , 
and (ii) there is apparently a small electron pocket at 
(0,0) 33 . Because our main intention is to understand 
what causes the pairing in the absence of hole pockets, we 
follow Ref. 26 and neglect the peculiarities of the folding 
procedure and potential electron pocket at (0,0). Cor- 
rections to this approach have been discussed in refer- 
ences 36 and 14. 

RPA/SF and fRG studies 26 ' 30 ' 31 applied to AFe 2 Se 2 
showed that the dominant instability is in the d-wave 
channel. The issue we want to address is whether the 
interaction between the two electron pockets alone is ca- 
pable to give rise to sizable d-wave attraction for large 
positive ^. We recall in this regard that for the two cases 
of electron doping that we considered before (/i = 0.05 
and fi = 0.18), the direct electron-electron interaction 
was attractive but very small, and the attraction in the 
d-wave channel was primarily the result of strong d-wave 
component of electron-hole interaction. Recall that at 
/i = 0.18 the hole FSs are already very small. Hole exci- 
tations are visible in ARPES in AFc 2 Se 2 above the gap 
of 60 — lOOmeV (Rcfs. 33), and it is certainly a possi- 
bility that the dominant mechanism of d-wave attraction 
between fermions on electron FSs in AFe 2 Se 2 is virtual 
pair hopping to gapped hole states. 

Let us elaborate on this point. When hole FSs are 
present, this pair hopping give rise to an effective at- 
tractive d-wave interaction between electron pockets with 
magnitude u^L, where L = logEp /T c comes from in- 
tegration over low-energy hole states. This interaction 
can well exceed a direct u ee because for relevant values 
of the logarithm Uh e L = 0(1). [A way to see this is 
to express A/, via A e in the matrix gap equation and 
re- write it solely as equation for A e .] When hole exci- 
tations are gapped, the logarithm is cut and scales as 
L = log \Ep/(Tc + Eq)] , where E is the gap for hole 
states. Because in AFc 2 Se 2 Eq < Ep (Ref. 33) it is not- 
apriori guaranteed that u e hL is now small and can be 
neglected. fRG study of superconductivity in AFe 2 Se 2 
(Ref. 30) does include this contribution because RG pro- 
cedure incorporates renormalizations coming from ener- 
gies above Ep. 

A way to verify whether the direct electron-electron in- 
teraction or the pair hopping to gapped hole states is the 
dominant mechanism of d-wave attraction in AFc 2 Se 2 is 



to use RPA/SF approach which only considers the in- 
teraction between fermions right at the FSs and neglects 
pairing interactions via intermediate gapped states, and 
see whether Ad is large enough (e.g., comparable to A s 
for hole doping), and whether the d-wave gap structure 
is similar to that obtained in fRG which includes virtual 
processes via gapped states. If Ad is not small and the 
gap structure is similar to that in fRG, electron- hole in- 
teraction is likely irrelevant, and d-wave attraction comes 
from the direct electron-electron interaction. 

In Fig. 12 we show the fit to electron-electron interac- 
tions obtained in RPA/SF formalism to Eq. 10 in which 
we only kept electron-electron interactions. The param- 
eters in s-wave and d-wave channels, extracted from the 
fit are shown in Table VII. 



TABLE VII: s- and d-wave parameters for the case of strong 
electron doping, when there are no hole FSs. 



s-wave 




OLee 


Pee X s 




0.84 


0.09 


0.04 -0.12 


d-wave 




OLee 


Pee Ad 




-0.04 


0.88 


-0.75 0.13 


TABLE VIII: s- 


and d-wave parameters for the case of strong 


electron doping, 


with band structure from Ref. 26, when there 


are no hole FSs. 








s-wave 




OLee 






3.65 


0.20 


0.03 0.1 


d-wave 




&ee 


Pee Ad 




-2.57 


0.29 


-0.0 5.9 



We solved 3x3 gap equations in s-wave and d-wave 
channels with these parameters, and show the results in 
Fig. 12. We clearly see that d-wave eigenvalue is positive. 

For completeness, we also computed the d-wave eigen- 
value for a different band structure, used in Ref. 26, which 
still has only electron FSs remaining. The results are 
shown in Fig. 13 and Table VIII. We find that the 
d-wave eigenvalue is again positive and the d-wave gap 
structure is quite similar to that in Fig. 12. The magni- 
tude of Ad, however, depends on the choice of the band 
structure. 

With LAHA fit, we are in position to analyze the pair- 
ing in more detail and check if the d-wave attraction 
comes from angle-independent or angle-dependent parts 
of electron-electron interaction. In the absence of hole 
FSs, Xd is the solution of 2 x 2 gap equation and is given 
by 

A d = -Uee (iPee + 1 ± ((2& e - l) 2 + 8a 2 £ ) ^ J (20) 

In general, Ad > either because u ee < 0, or u ee > 
but <5 2 e > Pee- In the first case, d-wave gap is due to a 
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FIG. 12: The fits of the interactions and the structure of s-wave and d-wave gaps in LAHA for the case of heavy electron 
doping, when only electron FSs are present. The As are X s — —0.12 and Ad = 0.13. (/i = 0.30) 
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FIG. 13: The same as in Fig. 12 but for a different band structure, used in Ref. 26. The couplings are A s =0.1 and Ad = 5.9. 



direct attraction between fermions from the two electron 
pockets, much in analogy with "hot spot" scenario for the 
cupratcs, and the gap has only modest variation along 
the electron FSs and no nodes. In the second case, d- 
wave gap should have nodes. We see from Table VII 
and VIII that u ee < , i.e. d-wave attraction in AFe2Se2 
is primarily due to the existence of a constant (angle- 



independent) attractive d— wave interaction between the 
two electron pockets. 

To understand why u ee < we remind that u ee 
is the difference between intra-pocket and inter-pocket 
electron-electron interactions u ee = (7Vp/2)(r eiei — 
r eie2 ) (see Eq. (10). The bare values of T eiei and 
r eie2 are quite close. For small electron doping, mag- 
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netic fluctuations are predominantly peaked at (0, ir) and 
(•7T, 0), and neither of these two interactions are selected 
for relative increase. As a result, the full intra-pocket 
and inter-pocket electron-electron interactions remain of 
nearly equal magnitude, and the difference between the 
two is small. Once the hole pockets disappear, the peak 
in magnetic susceptibility shifts to near (7r,7r), according 
to RPA 26 . In this situation, SF enhance the inter-pocket 
r eie2 compared to the intra-pocket r eiei , and u ee be- 
comes negative (attractive). 

We also note that s-wave solution remains attractive 
for one choice of the band structure (see Table VIII), 
although A s < A^. The two solutions for A s are given by 

A s = -u ee (2(3 ee + 1 ± ({2p ee - \f + 8o4) V2 ) • (21) 

The interaction u ee is given by u ee = (Np /2)(T eiei + 
r eie2 ) and is positive when both intra-pocket T eiei and 
inter-pocket T eie2 are positive, as in Tables VIII and 
VII). For u ee > a positive A s appears when a 2 ee > f3 ee 
which is satisfied for the parameters in Table VIII). 
Such a solution is induced by angle-dependent electron- 
electron interaction terms, and the s-wave gap should 
have nodes on the electron FSs. This is consistent with 
Fig. 13. 

The structure of both d— wave and s— wave gap func- 
tions is also quite consistent with the solution obtained 
within fRG 30 . The authors of that paper argued that the 
attraction in the d— wave channel is due to virtual hop- 
ping to gapped hole states, which is captured in fRG. We, 
on the contrary, believe that the good agreement between 
fRG and LAHA/RPA approach, which does not include 
gapped hole states, indicates that the primary cause for 
a d-wave pairing is the direct interaction between the 
two electron pockets, present in both approaches. We 
have not investigated the influence of gapped hole states 
directly, however. 

The issue that we don't discuss in this paper is how 
the d-wave gap evolves under the transformation into 
the folded zone. The folding is a non-trivial procedure 
for the case of d-wave pairing because the two electron 
FSs do hybridize in the folded zone, and this frustrates 
the d-wave gap which changes sign between the two un- 
hybridized FSs. Further complication is that AFe2As2 
has body-centered tetragonal lattice, and the two elec- 
tron FSs which eventually hybridize differ by k z = n and 
are rotated by ir/2 before hybridization. There is an ar- 
gument that in this situation the d-wave gap must have 
nodes near k z = 7r/2 (see Ref. 36). However, the issue 
of the pairing in the presence of the hybridization is not 
settled at the moment and we will not dwell on it. 

Finally, in our discussion we assumed that both intra- 
pocket interaction r eiei and inter-pocket interaction 
r ei e 2 are positive. Then u ee is definitely positive. Intra- 
pocket interaction is certainly positive, but the sign 
of interpocket interaction depends on the interplay be- 
tween U,V, J and J' in the same way as the sign 
of u\ le depends on the interplay between intra-orbital 



and inter-orbital interactions. If r eiC2 < the d-wave 
u ee = (Ap/2)(r eiei + |r eie2 |) becomes positive, while 
s— wave Uee = (NF/2)(T eiei — |r eie2 |) becomes negative, 
if |r eie2 | > r eiei . In this situation, the system develops 
an s— wave pairing with equal sign of the gap on the two 
hole pockets. We make a conjecture that r eie2 is nega- 
tive if one uses orbital J\ — J2 model with the interaction 
between first and second neighbors in real space (and 
J2 > Ji); instead of the on-site orbital Hubbard/Hund 
model. This would explain why both strong coupling 34,35 
and weak coupling 35 studies of the pairing in the orbital 
J\ — J2 model yielded an s— wave pairing. 



B. Strong hole doping 

Superconductivity at strong hole doping, when only 
hole FSs are present, has been observed in KFe 2 As2 40 ' 41 , 
which is at the opposite end from parent BaFe2As2in the 
family of Ki„, c Ba. E Fe2As2. T c in this material is rather 
low, only 3K, but the interest in KFe2As2is fueled by 
the fact that penetration depth and thermal conductiv- 
ity measurements point to nodal behavior, much as in 
LaFePO. The fRG study by Thomale et al. 29 shows that 
the leading instability is in d-wave channel, consistent 
with the observation of nodes. Interestingly, Thomale et 
al. have found that the d-wave gaps on the three hole FSs 
(two centered at (0,0) and one at (tt,tt)) are "in phase" 
with each other, i.e., A^(0) = Aj cos 2<fi (i = 1—3), and all 
Aj are of the same sign. This is not the d-wave solution 
with the largest A^ at smaller doping, when both hole and 
electron pockets are present. In that solution, there is a it 
phase shift between the gaps on the two FSs centered at 
(0,0) (see Figs. 10 and 9). The solution with "in phase" 
d-wave gaps is one of 5 d-wave solutions in the LAHA 
formalism, but with negative Xd at small/moderate hole 
doping. 

Given this discrepancy, it is interesting to analyze 
strong hole doping in our approach. We set fi = —0.3 
(n e = 4.88) and consider the same set of parameters as 
Thomale et al. (U = 3, J = 0.6, V = 1.85), but rescale 
all interactions by a = 0.25 to avoid a magnetic instabil- 
ity in the RPA/SF analysis. The magnetic susceptibility 
x(q) for such large negative \i has a maximum at q = 0, 
and affects mostly the interactions either near the FSs 
at (0,0) or at (jt,tt), but not intra-pocket interactions 
between (0,0) and (it, tt). 

The results for the fit of the RPA/SF interactions to 
LAHA arc shown in Fig. 14 and the interaction parame- 
ters extracted from the fit are presented in Table IX. We 
see that both s-and d-wave solutions are attractive and 
both are different from s— wave and d— wave solutions 
at smaller \/j,\, when hole and electron FSs were present. 
In particular, we see from Fig. 14 that s-wave gap now 
changes sign between the two hole FSs centered at (0, 0), 
while in the d-wave solution the gaps on these two FSs 
are "in phase". We remind the reader that at smaller \u\ 
the s-wave gaps on the two hole FSs are of the same sign, 
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while d-wave gaps on these FSs have a phase shift of tt 
(see Figs. 10 and 9). 

Consider first the s— wave channel. A simple analy- 
sis of the 3x3 equation shows that, in the absence of 
electron pockets, the reason for the s-wave attraction is 
a strong hi — h 2 interaction between the two pockets at 
(0,0), which exceeds the intra-pocket interaction for any 
of these two pockets (see Table IX). Intra and inter- 
pocket interactions are positive, and the solution with 
a positive A obviously corresponds to a sign-changing 
gap between the two hole pockets, in full analogy with 
the sign-changing solution in a model with strong inter- 
pocket interaction between hole and electron pockets. 
The presence of the third hole FS at (tt, tt) is not a factor 
in this consideration because the intra-pocket s-wave re- 
pulsion for this FS is quite strong, and the gap at (tt, tt) 
is only induced by much weaker interactions with (0, 0) 
pockets. 

Consider next the d-wave channel. A straightforward 
analysis of the 3x3 d-wave gap equation shows that the 
solution with A^ > exists for two reasons, both specific 
to the case of no electron pockets. First, we see from 
Table IX that the d-wave intra-pocket interaction within 
the (tt, tt) pocket, Uh 3 ,h 3 , is now negative (i.e., attractive). 
Second, intra-pocket interaction Uh t ,h 2 between the two 
pockets at (0, 0) is negative and larger in magnitude than 
repulsive Uh u h x and Uh 2 ,h 2 ■ In consequence, if we momen- 
tarily decouple the pocket at (tt, tt) and the two pockets 
at (0,0), we obtain two solutions with positive X d . One 
corresponds to a gap only on (tt, tt) pocket, another to 
in-phase gaps on the two pockets at (0,0). The third 
solution is the one in which there is the tt phase shift be- 
tween the two gaps at (0,0). This solution has negative 
Xd and is irrelevant. The residual, much weaker interac- 
tion between the pockets at (0, 0) and (tt, tt) couples the 
two solutions with positive Xd and sets the phase shift 
between the gaps at (tt, tt) and (0,0). In each of these 
two coupled solutions, there is no phase shift between the 
gaps at (0, 0), and the gap at (tt, tt) is larger than the gaps 
at (0,0) simply because Uh 3 ,h 3 , is attractive while Uhi,hi 
and Uh 2 ,h 2 are repulsive. These features are not present 
in the d-wave solution with the largest A^ for smaller 
when both hole and electron FSs are present. For those 
cases, Uh 3 ji 3 , is repulsive, and there is a tt phase shift 
between the two gaps at (0,0) because Uh x .h 2 is domi- 
nated by the interactions with electron pockets Uh ie and 
Uh 2 e which, by symmetry, are of different signs. We see 
therefore that the two d-wave solutions of the linearized 
gap equation which give the two largest Xd are not the 
same as the solution with the largest X d at smaller 

The larger value of A; l3 compared to A/^ and A^ 2 
and the in-phase structure of the gaps at (0, 0) are con- 
sistent with the fRG d-wave solution by Thomale et al. 29 . 
There is only one relatively minor disagreement: for our 
parameters the solution with the largest A^ is the one for 
which (tt, tt) gap and (0, 0) gaps have relative phase shift 
tt, i.e., are "of opposite sign" (sec Table X). Thomale et 
al. found the solution with the "equal sign" of all three 



d-wave gap. We verified, however, that the selection of 
the phase between (tt,tt) and (0,0) gaps is sensitive to 
the interplay between Uhi,h 3 and ilh 2 ,h 3 , which are small 
in magnitude and have different signs (see Table IX) . Al- 
ready a small modification of these parameters makes Xd 
larger for the solution with the same phase for the gaps 
on all three FSs, the same as in fRG solution. 

As about the comparative strength of s-wave and d- 
wavc pairing components, we found for the particular /x = 
—0.30 that we considered, that A s and Xd are comparable: 
A s = 0.13, X d = 0.11. However, the rate with which Xd 
increases with the hole doping well exceeds that for A s , 
and at larger dopings d— wave channel almost certainly 
becomes the most attractive one. This is also consistent 
with the fRG analysis by Thomale et al. who found that 
d-wave coupling becomes larger than s-wave coupling at 
large enough hole dopings. 

The pairing in heavily hole doped FcSCs was recently 
studied within RPA for 5-band orbital model by Suzuki 
et al 27 . They found that the pairing is driven by incom- 
mensurate spin fluctuations and s-wave and d-wave pair- 
ing amplitudes are of about the same strength. This fully 
agrees with our analysis. There is one difference, however 
- Suzuki et al attributed attraction in the s-wave chan- 
nel to the still strong interaction between hole states and 
gapped electron states near (0, tt) and (tt, 0), while in our 
case the attraction in s-wave channel is due to strong in- 
teraction between inner and outer hole pockets centered 
at (0,0) (see Table IX). The authors of Ref. 27 cited re- 
cent observation of incommensurate spin fluctuations in 
KFe2As2 (Ref. 42 ) as evidence for still strong interactions 
between fermions from near T and (tt, tt) points and from 
near (0,tt) and (tt, 0). This is certainly a possibility, but 
we point out that the interaction within the hole pocket 
centered at (tt, tt) also gives rise to incommensurate spin 
fluctuations at rather large momenta, because of a large 
size of that pocket. We recall that in our theory, the 
magnetically-enhanced interaction within the (tt, tt) hole 
pocket is the driving force for the d-wave pairing. 



C. strong electron vs strong hole doping 

We see that in our theory there is an attraction in d- 
wave channel at both strong hole doping and strong elec- 
tron doping (at least, for the model and parameters which 
we considered). The two limits are, however, quite dif- 
ferent from a physics perspective. For the case of strong 
electron doping, the enhancement of the spin susceptibil- 
ity around (tt, tt) unambiguously leads to a attraction in 
the d-wave channel, i.e. to Xd > 0. For strong hole dop- 
ing, the susceptibility is peaked at (0,0), which affects 
Uh lt hi, u hl ,h 2 , Uh 2 ,h 2 , and Uh 3 ,h 3 - The attractive d-wave 
solution is the result of negative Uh 3 ,h 3 and a larger value 
of h compared to u/i I ,/i 1 Uh, 2 ,/i 2 - There is no funda- 
mental reason why it should be so except that for large 
non-circular Fermi surfaces in 2D, the particle-hole sus- 
ceptibility x(l) is larger at 2k p than at q = 0. and the 
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appearance of the solution with a positive Ad at strong 
hole doping is very likely accidental. In other words, sim- 
ilar materials without electron pockets could easily have 
s-wave or a different d-wavc pairing state. 



TABLE IX: s- and d-wave parameters for the case of strong 
hole doping (with SF component), when there are no electron 
FSs. 



s-wave u hlhl 




Uh 3 h 3 


Uh 1 h 2 


Uhxh 3 


u h 2 h 3 


A s 


0.67 


0.69 


1.37 


0.80 


0.29 


0.48 


0.13 


d-wave u hl h 1 




Uh 3 h 3 


Uhxh.2 


Uh x h 3 




A d 


0.36 


0.34 


-0.17 


-0.50 


-0.02 


0.09 


0.11 



TABLE X: The structure of d-wave gaps Aj,^;) = 
Ah t cos 2<f>i for fi — —0.30 obtained by solving 3x3 linearized 
gap equation (the gaps are in arbitrary units since only the 
ratios of the gaps matter). 





sol. 1 


sol. 2 


sol. 3 




0.38 


0.60 


0.71 


A h2 


0.47 


0.54 


-0.70 


A h3 


-0.80 


0.60 


-0.08 


A 


0.11 


-0.06 


-0.43 



VI. CONCLUSIONS 

In this paper we analyzed the pairing symmetry and 
the structure of the gap in FeSC by approximating the 
pairing interaction between low-energy fcrmions by lead- 
ing angular harmonics. This allowed us to decompose 
the pairing interaction and study separately contribu- 
tions to pairing from scattering processes between dif- 
ferent FSs and the interplay between angle-independent 
and angle- dependent parts of each interaction. The an- 
gular dependence of the interactions are peculiar to FeSC 
because of the multiorbital nature of low-energy excita- 
tions. [The interactions in band representation are ob- 
tained by dressing up interactions in orbital represen- 
tation by angle-dependent coherence factors associated 
with the hybridization of Fe d-orbitals] . We used the 
band interaction obtained within RPA/SF formalism as 
an input, fitted it within LAHA, verified that the fit is 
quite good for all cases that we studied, and analyzed in 
detail how the pairing interactions in s-wave and c?-wave 
channels evolve with the bare interaction and the one 
with extra SF component and between hole and electron 
doping. We also analyzed the interplay between s-wave 
and d-wave pairing. Using the same procedure, we also 
studied the pairing at large electron (hole) doping, when 
only electron (hole) FSs are present. Throughout this 
paper we treated FeSC as quasi-2D systems and didn't 
address potential new physics associated with 3D effects. 



The main conclusion of our study is that all pairing 
states obtained so far at different dopings in FeSCs can 
be understood within the same universal pairing sce- 
nario based on spin-fluctuation exchange. We further- 
more found that all these pairing states appear naturally 
in the effective low-energy model with small number of 
input parameters. We conjecture that the approaches 
based on RPA (both analytical and functional) and on 
itinerant J\ — J2 model reduce to this model at low en- 
ergies, however with different input parameters. 

We used this effective model to study the doping evo- 
lution of the pairing in hole and electron-doped FeSCs. 
We argue that the pairing mechanisms at small/moderate 
and large dopings are qualitatively different - when both 
hole and electron pockets are present, the pairing is 
of Kohn-Luttinger type, driven by the pair-hopping of 
fcrmions from hole to electron pockets, while at larger 
hole or electron doping, the pairing is due to a direct in- 
teraction between only hole or only electron pockets. For 
moderate hole dopings the leading pairing instability is 
towards an s± state with a nodelcss gap. For moderate 
electron doping a nodal s± state is the leading instabil- 
ity, but d-wave is a close competitor. For larger electron 
or hole dopings, when only one type of FS is present, 
the leading pairing instability is towards a d-wave state, 
which in case of strong electron doping is nodeless, at 
least in 2D case. 

We summarize below the detailed reasoning behind the 
observation stated above by presenting the answers to the 
questions we posed in the introduction- 

• What is the origin of the strong angular dependence 
of the s^ gap along the electron FSs? 

We found that the origin is different for bare and 
full interactions. For bare interactions (no SF 
component), the combination of intra- and inter- 
electron pocket repulsions are stronger than the 
electron-hole interaction and the s^ attractive so- 
lution for the gap is entirely due to the angle- 
dependent parts of the electron-hole (electron- 
electron for strong electron doped materials) inter- 
action, much as was anticipated in Rcfs. 1 and 
43. Namely, the system adjusts the magnitude of 
the angle-dependent, ± cos 28 gap component along 
the two electron FSs to minimize the effect of the 
inter-electron-pocket repulsion. 

For the full interaction, the electron-hole interac- 
tion is the strongest, and the attractive s^ so- 
lution exists even if all interactions are taken to 
be angle-independent. The angle-dependent terms 
modify the s^ gap by creating ± cos 29 gap compo- 
nents. Whether these components are large enough 
to lead to nodes depends on details, but the generic 
trend is that when the angle-independent part of 
the electron- hole interaction is larger, the gap is 
less likely to have nodes. 

• Are the angular dependencies of all interactions rel- 
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FIG. 14: The fits of the interactions by LAHA and the structure of s-wave and d-wave gaps for the case of heavy hole doping, 
/i = —0.30, when only hole FSs are present. The parameters are U = 0.75, J = 0.15, V = 0.46 (see text). The eigenvalues are 
\ a = 0.13 and \ d = 0.11. 



evant for the gap structure, or can some interac- 
tions be safely approximated as angle-independent? 

We found that the angle-dependent part of the 
electron- hole interaction is the relevant one. The 
angle-dependent parts of the electron-electron in- 
teraction have little effect on the gap structure, at 
least for the full interaction with the SF compo- 
nent. 

• Why do the solutions obtained within the 
RPA/SF and fRG approaches have nodes for sys- 
tems with two hole and two electron FSs and no 
nodes for systems with three hole and two electron 
FSs? 

We found that the angle-independent electron-hole 
interaction, which favors a no-nodal s^ gap, is fur- 
ther increased if the third hole FS is present. For 
most of the parameter sets which we analyzed, the 
gap has nodes in case of electron doping (four 
FSs), but no- nodal solution is stabilized for hole 
doping (five FSs). Kuroki et al 3 have pointed out 
that this can be traced back to the d xy orbital 
character of the third hole pocket which interacts 
strongly with the d xy states at the tips of the elec- 
tron pockets. We found, however, that for some 
parameters nodal solutions survive in the presence 
of the fifth FS, i.e., the disappearing of the nodes 
with the appearance of the fifth FS is not a uni- 
versal result. That aside, the gap structure still 
evolves between nodal and no-nodal once we change 
the magnitudes of angle-dependent parts of the in- 



teractions. 

• What causes the pairing when only electron FSs are 
present? 

We found that the d-wave pairing is generally at- 
tractive and competes with s^ pairing for the 
electron-doped FeSC. At small electron doping, the 
(i-wave attraction is almost entirely due to the d- 
wave electron-hole interaction, and the direct d- 
wave interaction between electron pockets is weak. 

For strong electron doping, when only electron FSs 
are present, the situation is different. We found 
an attractive d-w&ve interaction between electron 
pockets. The d-wave pairing is then quite similar 
to the one in the magnetic hot-spot pairing sce- 
nario for the cuprates. In both cases, there is a 
(i-wave attraction between the FS sheets separated 

by (7T,7r). 

With regard to the sublcading s-wave attraction, 
in our case it is due to an angle-dependent s- 
wave component of the electron-electron interac- 
tion. We didn't consider the interaction via gapped 
hole states (another potential reason for s-wave at- 
traction), but the similarity between our gap struc- 
ture and the one obtained in the fRG study 30 , 
which includes both interactions, indicates that the 
likely origin of the s-wave attraction is the angle- 
dependence of electron-electron interaction. 

• What causes the pairing when only hole FSs are 
present? 
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We found that both s-wave and d-wave channels 
are attractive, with comparable X sd > 0. Which 
pairing instability is stronger depends on detail. 

The reason for the s-wave instability is a strong re- 
pulsive inter-pocket (hi — h 2 ) interaction between 
the two pockets at (0,0), which exceeds the intra- 
pocket repulsion. This leads to a sign-changing 
s-wave gap between the two hole pockets, in full 
analogy with the sign-changing solution in a model 
with strong inter-pocket interaction between hole 
and electron pockets. 

The reason for the competing d-wave instability is 
two-fold. First, the d-wave component of the intra- 
pocket interaction within the hole pocket at (%, tt) 
is negative (i.e., attractive), second, there is strong 
attractive d-wave interaction between the two hole 
pockets at (0,0). The combination of these two 
reasons leads to a d x i_ y i solution with positive Xd, 
in which the magnitude of the gap is the largest on 
(tt, tt) pocket, and the two gaps at (0, 0) have zero 
phase shift. 

Both s-wave and d-wave solutions are different from 
the ones at smaller hole dopings, when hole and 
electron FSs are present (e.g., the cZ-wave solution 
with the largest Xd at smaller \fj,\ is the one with 
the tt phase shift between the two gaps at (0,0). 

• How is the structure of the pairing interaction af- 
fected when the spin-fluctuation component is added 
to the direct fermion-fermion interaction? 

We found that the SF interaction primarily changes 
the overall magnitude of the interaction, while 
its angular dependence remains nearly unchanged. 
All components of the pairing interaction increase 
when the SF term is added. On top of this, there 



is an additional increase of the hole-electron inter- 
pocket interaction, both in the s-wave and in the 
(i-wave channels. This additional increase makes 
both s-wave and d-wave solutions attractive. 

We have only studied the strictly 2D case thus far, and 
neglected aspects of the 3D I4/mmm crystal symmetry 
characteristic of 122 materials and the hybridization of 
electron pockets in the folded zone. We nevertheless be- 
lieve that the general evolution of interactions and gap 
symmetry discussed here will be generic to the FeSCs. 
The approach developed here can be easily modified to 
study superconductivity on hybridized FSs and can also 
be used to study in great detail SDW instability in multi- 
orbital systems 44 . 
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